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We study the presence of chaotic behavior in phase space in the pre-inflationary stage of hybrid 
inflation models. This is closely related to the problem of initial conditions associated to these 
inflationary type of models. We then show how an expected dissipative dynamics of fields just before 
the onset of inflation can solve or ease considerably the problem of initial conditions, driving naturally 
the system towards inflation. The chaotic behavior of the corresponding dynamical system is studied 
by the computation of the fractal dimension of the boundary, in phase space, separating inflationary 
from non-inflationary trajectories. The fractal dimension for this boundary is determined as a 
function of the dissipation coefficients appearing in the effective equations of motion for the fields. 
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. I. INTRODUCTION 

> . 

Among the various proposed models for implementing inflation, the hybrid inflation model |l| is one of the most 
' attractive ones due to the possibility of its implementation in the context of supersymmetry and supergravity models 
|§| • The principle behind these models for inflation are basically based on models where the inflaton is coupled to one 
' or more scalar fields. The inflationary phase is characterized by an initial phase in the evolution of the fields where 
T-H , the inflaton field slowly moves towards zero vacuum expectation value, till it approaches a critical value which then 
induces a spontaneous symmetry breaking in one of the other fields which are coupled to the inflaton, after which the 
fields quickly evolve to their vacuum states and infiation ends. 
^ I It has been shown in the context of supergravity motivated models that for inflaton field amplitudes larger 
' than the Planck mass there will appear large quantum corrections in the inflaton's effective potential, destroying the 
^ . flatness of the inflaton's potential as required by the slow-roll conditions for inflation. The question then turns to 
' whether inflation can be achieved in these models for initial inflaton field amplitudes smaller than the Planck scale, 
. where a reliable model could be constructed. But in this case it has been shown |^ that small fluctuations of the other 
^ ' fields coupled to the infiaton are efficiently able to prevent the onset of infiation, or quickly drive previous infiationary 
• I— I , trajectories to non-inflationary trajectories and therefore making inflation to end before enough e-folds of inflation 
' have been produced [Ng ^ 70 as required) to solve the usual cosmological horizon and flatness problems. This is 
\ related to the homogeneity requirement over Hubble size distances ||^ in order to inflation to begin. In the context of 
- - hybrid inflation this problem is even more severe, requiring an extreme fine-tuning of the initial conditions in order 
to have sufficient homogeneity, with negligible spatial and time derivatives of the other fields that are coupled to the 
infiaton, over regions extending far beyond a few Hubble lengths. 

A few solutions to this homogeneity, or fine-tuning problem, have been proposed. The authors in Ref. Q for 
example propose two-stage infiation models where a first short inflationary phase would smooth a large enough region 
in phase space making then possible a later longer second inflationary phase. In Ref. it is analyzed the possibility 
of the solution of the fine-tuning problem in the context of more complicate generalizations of the hybrid infiation 
model as motivated from higher-dimensional scalar field models or brane cosmology. 

Here I will be interested in investigating the possible role of how effective field interactions may shed further light 
in this fine-tuning problem in the pre-infiationary phase. This is mostly motivated from several works on the effective 
dynamics of scalar field configurations, which have shown that dissipation is intrinsically present in the field dynamics 
|p|-|l5|. In special, the role of dissipation in all of these works has been emphasized. We can trace the origin of these 
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effects by considering the dynamics of a given system, which interacts with a sufficiently large environment (in the 
sense of degrees of freedom). One well known example of this, in the context of statistical mechanics, is the model of 
an oscillator, taken as being the system, in interaction with a large number of harmonic oscillators which are taken 
as representing the environment. By functionally integrating over the environment degrees of freedom we can then 
show that the dynamics of the system oscillator will be dissipative, with an equation of motion of the form of a 
Langevin-like one 

The dynamic of interactingfields has being a subject of intense study, in both Minkowski space-time [p| Jl^Jl5[ |, as 
also in FRW space-times |]lT| , ]l2| , with special interest in studying the dynamics of the inflaton field both during and 
immediately after the infiationary phase. In particular, in Refs. [ p"7[p^ it was shown the importance of considering 
dissipative effects (due to inflaton's decay) during inflation. The authors in have shown how dissipation can 
influence the usual scenario of inflation, as also the effect of dissipation on the field trajectories in phase space. Strong 
dissipative regimes for field evolution have in special motivated the implementation of new inflationary models called 
"warm inflation" . The viability of the construction of microscopic motivated models of strong dissipative inflaton 
models has been shown possible in the context of superstring inspired models in Ref. |^ and argued also possible for 
more general models in [|2^ . 

On very general grounds, we are then lead to inquire about the possible role of field dissipation also at the early 
stages of inflation and in special to its role in smoothing large regions in phase space and thus easing considerable the 
homogeneity requirements for the onset of inflation and, consequently, providing a natural solution for the problem 
of fine-tuning associated to hybrid inflation. It is easy to understand the effect of dissipation on the initial conditions 
problem. For the regime of initial inflaton amplitude below the Planck scale the energy density is small at the 
beginning of inflation, resulting in a small Hubble parameter H , which determines how fast the inflaton energy is 
converted to expansion. For small values of H the coupling of the fields to the background metric of gravity is small 
and the corresponding friction like term in the equations of motion, of the form 3H(f), is small. This results in a very 
long evolution for the fields, where those fields coupled to the inflaton oscillate around zero several times, generating 
strong sensitivity to very small variations of the initial conditions and resulting in a complete indetermination of the 
flnal outcome of the fields' trajectories, which can be evolution toward an infiationary regime or to the minima of 
the potential. The behavior of the fields' trajectories in the early times before a required long inflationary regime 
is then extremally chaotic. In the language of dynamical systems we have two types of attractors in the system, 
represented by the inflationary regime and the minima of the potential. The presence of dissipation works in damping 
the fluctuations of the fields in this initial critical time, suppressing the chaotic motions, which we then expect, from 
the results of Ref. [0, will bias the inflaton field trajectory toward the inflationary one. In fact, recently a study done 
by the authors in Ref. [ p2| have indicated that additional damping terms in the inflaton equation of motion could 
alleviate many of the problems related to the homogeneity requirements before inflation. Another study done by the 
authors in Ref. [23| , using a supersymmetric hybrid inflation model, have also indicated that particle production may 
be a way to relax the extreme homogeneity requirements for hybrid inflation. 

Here we will then be mostly concerned in determining for what magnitude of dissipation will the outcome of the 
evolution of the fields tend mostly toward the inflationary regime, or in the dynamical system approach, determine 
how fast dissipation will change the chaotic regime, characterized by unstable inflationary trajectories, to a non- 
chaotic one, with stable inflationary trajectories. We study chaos in our dynamical system of equations by means 
of the measure of the fractal dimension (or dimension information) [for a review and deflnitions, see e.g., Ref. p^], 
which gives a topological measure of chaos for different space-time settings and it is a quantity invariant under 
coordinate transformations, providing then an unambiguous signal for chaos in cosmology and general relativity 
problems in general ||25| , ^6|| . The method we apply in this work for quantifying chaos is then particularly useful in this 
cosmological pre-infiationary scenario context we are studying, in which case other methods may be ambiguous, like, 
for example, the determination of Lyapunov exponents, which does not give a coordinate invariant measure for chaos, 
as discussed in [ 25|j2^ . Also, other methods for studying chaotic systems, like for example by Poincare sections, are 
not suitable in the case we are interested here, in which case chaos is mostly a transitory phenomenon (it ends by the 
time the fields reach the potential minima, or when the infiaton enters the inflationary region). By computing the 
fractal dimension characterizing this chaotic behavior, which is related to the uncertainty in the system parameter 
values to predict the final outcome from a given initial condition (we may say that the boundaries of initial conditions 
that lead to infiation or evolution towards the potential minima are mixed) we are then able to infer the naturalness 
of inflation for different settings of the initial conditions. To our knowledge, this is the first time that this kind of 
study is performed in the context of a dissipative dynamical system in cosmology and applied to the initial condition 
problem of inflation in particular. 

The paper is organized as follows: in Sec. II I give the basic equations deflning the dynamical system and I discuss 
some of their properties. In Sec. Ill I present the numerical analysis of the dynamical system and the computation of 
the fractal dimension as a function of the magnitude of dissipation. From this analysis one will be able to conclude 
about the general effect of dissipation in the evolution of the fields and how it works in favor of the inflationary regime. 
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Finally in Sec. IV I give the concluding remarks. 



II. THE MODEL AND ITS PROPERTIES 



The model I will study here consists of the simplest hybrid inflation model, with a scalar field (f> (the inflaton) 
coupled to another scalar field cr, which triggers the end of inflation. The potential is 

1/(0, a) = jia'- M^f + 1^02 + , (2.1) 

where the parameters M^, rn? and the couplings A and arc positive. In the following I will treat only the case of 
homogeneous fields, = (j){t) and a = cr(t), which is the usual case in the works dealing with the initial condition 



problem in hybrid inflation. The interpretation of inflation from the potential (2.1) is the standard one. For values 
of (j) larger than a critical value 0cr, where 0^j. = \M'^ / , there is no symmetry breaking in the cr-field direction and 
tr = is a local minimum of the potential. We are interested in the region where Mpj ^ 4>^ ^ 0cr where inflation 
takes place (the false vacuum dominated regime). After the inflaton field drops below 0cr, symmetry breaking occurs 
in the a field direction of the potential. At this point the fields quickly move towards the minima a = ±Af , = 
and inflation ends. 

The chaotic properties of the classical (homog eneo us) equations of motion, in Minkowski spacetime, for a model 



with potential similar to the one given by Eq. (2.1) were studied in Ref. JB7|, while the full effective equations of 
motion in Minkowski spacetime were studied by Ramos and Navarro in Ref. p8|]. In Ref. |2^ the general form expected 
for the equations of motion for the fields, for potentials of the form of Eq. ( |2.1D , was derived and a detailed account 
for the dissipative terms appearing in those equations was given. The chaotic behavior of the equations of motion for 
the fields was quantified by means of the fractal dimension and the authors studied in details how dissipation of the 
fields, due to decaying modes, changed the chaotic behavior of the dynamics. I will here then extend the results of 
Ref. psl to the case of an expanding background, a flat Friedman-Robertson- Walker background metric. 



As shown in Refs. [^j28|2l[] dissipation comes from the coupling of the system fields and a (here taken as 
background field configurations) to a bath made of a set of other fields, for example made of fermions ipi and/or 
N-)^ scalars Xi- The general form of the interactions can be of the form 

T^ial^'xl + gWxj) , (2.2) 

i=l 

for scalar fields Xi and 

^{f<t,4'i(l)ipi + fa-ipicril^i) , (2.3) 

i=l 

for fermion fields tpi- The usual form expected for the coupled effective equations of motion for and cr are complicated 
non-local equations of motion with typical non-Markovian dissipative kernels. In ||2^ we have shown that at high 
temperature and in the large N-^, limit for the (thermal) bath fields we can find an approximate Markovian 
limit for the kernels, from which we can express the equations of motion for the fields in terms of dissipative local 



equations. In Ref. |21| we have studied in details the general equations for the non-Markovian dissipative kernels at 



zero temperature and shown that for a certain class of field decaying modes the Markovian approximation is also 



a valid assumption. The results obtained in Ref. |21| show that we can evade the assumption of an initial high 
temperature thermal bath as used in Ref. to derive the local equations of motion for the fields and makes it, 
therefore, more appropriate to the application we intend for here, which is the study of t he in itial conditions problem 



for the onset of inflation in hybrid inflation models characterized by potentials like Eq. (2.1). 

I will just use the general form of the coupled effective equations of motion obtained in Ref. without a complete 
specification of the dissipative coefficients, which we will take as free parameters. This is a valid analysis here since 
I will mostly be interested in how the general dissipation expected for the fields and a in their coupled effective 
equations of motion will change the chaotic properties of the dynamical system, which, as explained in the introduction, 
are direct related to the initial conditions problem. The detailed form of the dissipative coefficients depends on the 
microscopic physics of the specific model under study, like the coupling of the system fields with the bath fields, the 
decaying modes available and the expanding metric. In fact the magnitude of the dissipation terms may be controlled 
by these couphngs, as shown in Refs. [ p^pl|j2^ . 
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In analogy to the results obtained in Ref. |28) we write the equations of motion for and a in the form|^ 



(2.4) 



and 



(2.5) 



where rji, 772 and 773 denote the dissipation coefficients. Besides the coupled equations of motion (2^) and ( |2.5| ) we 
also have the Friedman equation and the evolution equation for the Hubble parameter, H = a/a, where a is the scale 
factor, are given as usual by 



SttG 



[Pm + Pr) 



(2.6) 



(2.7) 



where G — 1/mpj, with mpi the Planck mass, k = 0, +1, —1 for a flat, closed or open Universe, respectively. Pm{r) ^nd 
Pm{r) are the energy density and pressure for matter (radiation), respectively. We also have the standard relations: 



(2.8) 



and Pr — ^Pr- 

The matter and radiation energy densities Pm and Pr evolve in time as: 



(2.9) 



(2.10) 



and (from the energy conservation law) 



Pr + 4:Hpr — rii(f)^(p — ri2<y^cr^ — 2i]:j(f>a4>& ~ 



(2.11) 



Assuming a flat universe (fc = 0), from Eq. (2.6), we can consider 



Pr = 



8ttG 



8^^ -Y-Y-^(^'^) 



(2.12) 



as the first integral of Eq. (2.11). Using Eqs. (2.6)-(2.9), we can also express the equation for the acceleration in the 
following form 



a _ SttG 
a 3 



(2.13) 



The field equations (2.4) and (2.5) with ( ^.13 ) form a dissipative dynamical system that we will study numerically. 
Using ( 2.12| ) in (2.13) and defining the following dimensionless variables: y = VSvrGi/), x — cr/M, r = Mt and the 
constants = rnTjM'^, 0^ — SttGAP an d the rescaled dissipative coefhcients rjy — Mrji, rjx — Mri2 and rj^y — Mrj^, 
the system of equations (^.4|), (2.5) and ( 2.13 ) can then be rewritten in terms of dimensionless variables in the form 
of the following system of first order differential equations: 



^This should be compared with Eqs. (5.1) and (5.2) of Ref. [E8|. In the notation used here we have that in those equations 



that 0c 0, -ipc 



AM , — > 0, — > 6A and g ^ g. 
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y' = w , 
a = u , 

z' = -3^z + \x ~ ~ §2 ^2/^ - rjxX^z - -^'HxyXyw , 
w' = -3^w - a'^y - g^yx^ - ■^r/yy'^w - ijxyxyz , 

u' = - - (/J^z^ + _ _^ 2A/32x2 - A/32a;4 - 2a^y^ - 2g^x^y^) , (2.14) 
a 6 

where prime indicates derivative with respect to the dimensionless time r, e.g., x' = dx/dr. 

For convenience we use the values for A, g^, M and m as given by Mendes and Liddle in Ref. 0: X = g^ = 1, 
AI — 2 X 10~^Mpi and m = 5 x lO^^A/pi, where M pi = I/VSttG is the reduced Planck mass. From these values of 
TO and M we find for the constants a and /3 in Eq. ( ^.14[ ) the values: a — 2.5 x 10^'' and f3 = 2 x 10~^. In terms of 
these values we also have that the critical value for the inflaton field, (/)cr, in the dimensionless variables, is given by 
ycr = 2 X 10-2. 



III. CHAOS AND THE FRACTAL DIMENSION 



We next numerically solve the system of equations in Eq. (2.14) and we search for chaotic regimes and how they 
change with increasing dissipation. Lets initially consider the following initial conditions at i = 0: = 4 x (per, = 0, 
cr = 2 X IQ-^M, & — and H is determined initially by Eq. ( |2.6| ), with zero initial radiation energy density. These 
values guarantee that the stable inflationary trajectories will correspond to at least ^ 60 e-folds of inflation. For 
the dissipation coefficients rj^, rjy and rj^y, for convenience, we will consider them all with the same value, which is 



consistent with the recent findings in Ref. |2lf| , or the calculations in Ref. 28 1 . Small variations of the magnitude among 



these coefficients are not critical here, since the dissipation in Eq. ( 2.14 ) is dominated by the term corresponding to 
rji(j>^(f>. We take the dissipation coefficients with initial values as given by {77} — {rjx , ^ Vxy) = (0,0,0) and increase 
them till the chaotic behavior of the system changes to non-chaotic. 

For {rj} = a typical result for three trajectories around the initial values of the fields (in dimensionless units), 
separated by a ±10"^ variation, is shown in Fig. 1. 

In Fig. 1 we can distinguish the inflationary trajectory represented by the straight line at x — 0, which represents 
the evolution of the fields in phase space along the valley of local minima of the potential, cr = for > (per. This 
trajectory is characterized by a very long evolution by which the universe expands exponentially. The number of 
e-folds produced by the end of the inflationary phase, when the inflaton reaches (p ^ </'cr, for this particular trajectory 
is Ne ^ 350. The other two oscillatory trajectories, that evolve towards x — ±1 represent the non-inflationary 
trajectories, where the fields quickly evolve towards the potential minima a = ±M and (f) — 0. In Fig. 2 we show a 
blow-out of the initial time evolution of Fig. 1. It clearly shows highly oscillatory chaotic like behavior. 

The chaotic behavior of the dynamical system is quantified by means of the determination of the fractal dimension 
of the boundary separating the inflationary trajectories from those that are not infiationary, i.e., the ones that evolve 
towards the minima of the potential, (p — and a = ±M. The fractal dimension is associated with the possible 
different exit modes under small changes of the initial conditions at < = and it will give a measure of the degree of 
chaos of our dynamical system. The exit modes we refer to above are ones of the symmetry breaking minima in the 
tr-field direction, a = ±M, (p = and the inflationary one, which are attractors of field trajectories in phase space. 
The method we employ to determine the fractal dimension is the box-counting method which is a standard method 
for determining the fractal dimension of boundaries [2^ ] . Its definition and the specific numerical implementation we 
use here have been described in details in Ref. |^ . 

The basic procedure is, given a set of initial conditions xq at i = ioi which leads to a certain outcome for the 
trajectory in phase space, by perturbing them by an amount i5 we then study whether there will be a change of 
outcome for the trajectory or not (whether the perturbation will lead to a different attractor or not). Given a volume 
region in phase space around a boundary between different attractors and perturbing a large set of initial conditions 
inside that region, the fraction of uncertain trajectories, f{S), which result in a different outcome under a small 
perturbation, can be shown to scale with the perturbation S as f{5) ^ S'^, where e is called the uncertainty 
exponent. The box-counting dimension of the boundary in phase space separating different attractors, or fractal 
dimension fd, is given by ||24| fd — d — e, where d denotes the dimension of the phase space. For a fractal boundary 
fd > d — 1, implying that e < 1, whereas for a non- fractal boundary, fd = d — 1, and e = 1. 
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Following the method of box-counting we then consider a box in phase space (for the dimensionless variables) of 
size 10""', around the initial conditions used in Figs. 1 and 2, inside which a large number of random points are taken 
(a total of 100.000 random points were used in each run). All initial conditions are then numerically evolved by using 
an eighth-order Runge-Kutta integration method and the fractal dimension is obtained by statistically studying the 
outcome of each initial condition for each run of the large set of points. Special care is taken to keep the statistical 
error in the results always below ~ 1%. From these numerical simulations we obtain for the zero dissipative regime 
of Fig. 1 the result for the fractal dimension as given by fd — 5.80 ± 0.05. This corresponds to the dimension of 
the dynamical system, Eq. (2.14), corresponding here to d = 6, minus the uncertainty coefficient e (which gives a 
measure of how chaotic is the system [^), /d = 6 — e. The results obtained by increasing the value of the dissipation 
coefficients {ry}, while keeping the same initial condition used above, around which the perturbations are taken, are 
shown in Fig. 3. 

The results in Fig. 3 show that the system quickly changes to a non-chaotic behavior for a dissipation coefficient 
around {77} ~ 3.1 x 10"''. For those values of dissipation and higher all trajectories are inflationary ones with ^ 60. 
At that point we then have the breakdown of the fractal structure of the boundary between inflationary and non- 
inflationary trajectories. The boundary becomes smooth and we can reliable predict that the initial conditions will 
evolve towards the inflationary region. 

In order to study how efficient is dissipation to drive non-inflationary trajectories towards the inflationary region, 
which is characterized by a > 0, we show in Fig. 4 the results for the number of e- folds of inflation as a function of 
dissipation for an initial trajectory that in the absence of dissipation would be a non-inflationary one. The particular 
initial conditions we take purposefully correspond to fields well below the Planck scale, 4> — 3.5M, cr — M, and 
non- vanishing field derivatives (p = 2.5 x 10~^M^ and o- = 5 x 10~''M^, which in the absence of dissipation would be 
quickly driven to one of the minima of the potential (and therefore corresponding to a non-infiationary behavior). We 
see that above some value of dissipation coefficients 77 ~ 1 the number of e-folds of inflation increases fast for rather 
very small increase of dissipation. We then clearly see using the set of initial conditions above how dissipation acts 
in turning non-inflationary trajetories into inflationary ones, shifting the system deep into the inflationary region as 
dissipation is increased. The initial radiation energy density in this case, as in the previous one, is taken initially zero 
and remains always smaller than the vacuum energy density during the whole evolution of the inflaton. In particular, 
for the values of dissipation coefficients we have studied here, after a few e-folds of inflation the radiation energy 
density is a negligible fraction of the vacuum energy density and right before the end of inflation the radiation energy 
density is a tiny fraction of the vacuum energy density, Pr/Pv ^ 10"^**. A detailed treatment of the liberated radiation 
to the medium and its eventual translation to a temperature depends of course on how the radiation thermalization 
rate Fj-ad, which is determined by the microscopic physics, compares to the expansion rate H. 

It is also useful to compare the regimes of field dissipation corresponding for instance to the ones studied in Fig. 
3 and Fig. 4 and their relation to the warm infiation scenario. In warm inflation ||r^ , |20|| we are usually interested in 
the regime of strong dissipation in the effective equation of motion for the inflaton, as compared with the expansion 
term 3H(j). From Eq. (|]J) dissipation is dominated by the term rii(jp(j) ~ ri[cj))(j). From the values of parameters 
and initial conditions used in Fig. 3 we flnd that 3H/r]{(f)) 5 for the largest value of dissipation studied, while 
for the parameters and initial conditions used in Fig. 4 we have that ZH/rj{(f>) ~ 0.008 and therefore deep in the 
regime of warm inflation. We can easily see then that warm inflation kind of models does not seem to suffer from 
the initial conditions problem seen in standard hybrid inflation models. In particular, for the parameters used in Fig. 
3, for values of dissipation 77 > 0.0015 we enter in the warm inflation regime and the system has well before became 
non-chaotic, with stable inflationary trajectories. 



IV. CONCLUSIONS 



We have examined a dynamical system describing hybrid inflation in the presence of dissipation for both the inflaton 
and the field that triggers the end of inflation. Dissipation is taken as being present throughout the system's evolution. 
In the absence of dissipation the dynamical system is shown to be highly chaotic by means of the measure of the fractal 
dimension. The fractal dimension gives a coordinate invariant measure for chaos and therefore is an appropriate way 
to quantify chaos in cosmology. 

For a given initial condition that evolves to an inflationary region, by varying it by a small perturbation we 
completely loose track of its final outcome, which may still evolve towards the inflationary region, characterized by 
evolution along the valley of local minima of the potential, characterizing the inflationary regions, or fast evolution 
towards one of the symmetry breaking minima. The time evolution of the fields are shown to be highly oscillatory 
for a very long time, which eventually change the inflationary trajectories to non-inflationary ones. This is the main 
characteristics in the evolution of the fields and that gives rise to the initial condition problem in hybrid inflation 
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(and in a lesser degree in any inflationary model). The long oscillatory behavior of the fields is highly chaotic. By 
measuring the magnitude of chaos in the system we can then infer the severity of the initial conditions problem. 

However, by changing the magnitude of dissipation we have then shown that the highly chaotic evolution of the 
system, with unstable inflationary trajectories, can change quickly to a non-chaotic, stable one, where all trajectories 
eventually evolve toward the inflationary region. Dissipation acts by damping the initial oscillatory behavior of the 
fields and then making it possible to provide a solution for the initial condition problem. We have shown that this 
solution can be obtained with relative very small dissipation of the initial inflaton amplitude, with rates of conversion 
of vacuum energy to radiation as small as ~ 10^^*. Here dissipation works in two ways. By damping the fluctuations 
associated to the fields coupled to the inflaton the fine-tuning problem associated with hybrid inflation models is 
avoided and at the same time it makes possible for initial inflationary trajectories to enter faster in the inflationary 
region (characterized by a > 0) and stay longer in there. 

The method used here to study the initial condition problem can also be easily extented to include other variants of 
the hybrid inflation model, for example in the context of supersymmetric motivated potentials, where more than one 
field is coupled to the inflaton. The additional equations of motions will also generically be dissipative equations and 
the same changing of behavior due to dissipation observed here is also expected to appear in those more complicated 
models. 

There have been previous studies on chaos in hybrid inflation models ^o|-^ , but these studies have concentrate on 
the reheating period after the inflationary phase. In that context, the effect of dissipation, that would be inherently 
present in the field equations also during this period, has been neglected. It would then be interesting to investigate 
in those cases also the effect of dissipation and its role, for example, in the phenomenon of parametric resonance 
during preheating in hybrid inflation models p^ , |3^ . We are currently investigating this and other applications and 
we expect to report on them soon. 
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FIG. 1. A typical evolution of three trajectories in the absence of dissipation with initial conditions separated by a ±10 ^ 
variation around the t — values </> = 4 x ipcr, (j> = 0, a — 2 x 10~^M, a = and H determined by Eq. (2.6). 
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FIG. 2. A blow-out of the initial time evolution of the trajectories shown in Fig. 1. 
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FIG. 3. The fractal dimension as a function of dissipation, from the box-counting method appUed to a box in phase space 
centered around the initial conditions used in Fig. 1. 



300 



250 



200 




150 



100 



0.91724 



0.91726 



0.91728 



0.91730 



FIG. 4. The total number of e-folds of inflation as a function of dissipation for the initial conditions at t = 0: 
(T = M, = 2.5 X 10-^M^ CT = 5 X IQ-^M^ and pr{t = 0) = 0. 
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